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Abstract 

We extend some results of Itai Benjamin! and Yuri Lima (see Tl). In this paper they 
consider a binary tree T„ of height n, each leaf is either infected by one of k diseases or 
not infected at all. In other words, x at generation n is infected by the i-th infection with 
probability pi and sane with Pk+i- Moreover the infections are independently distributed for 
each leaf. 

Infections spread along the tree based on specific rules. In their paper they study the limit 
distribution of the root of T„ as n goes to infinity. 

Here we want to study the more general case of a Galton-Watson tree and a 2-ary tree. 


1 Introduction 

First we recall the definition of a Galton Watson tree (GW) and give a few notations. Assume 
that is a N-valued random variable following a distribution q, in other words ¥{N = i) = qi for 
f G N and to have an interesting problem, we assume that go + <Zi = 0 (Botcher case). 

Let (j) be the root of the tree and an independent copy of N. Then, we draw children of 
(j): these individuals are the first generation. In the following we write N for for typographical 
simplicity. At the n-th generation, for each individual x we pick an independent copy of N 
where is the number of children of x and so on. The set T, consisting of the root and its 
descendants, forms a GW of offspring distribution q. 

We denote by |a;| the generation of x and for n S N, T„ = {x G T, |a;| < n} the GW cut at height 
n and the leaves of T„ are the elements of T„\T„_i. 

According to Neveu’s notation ([3]), to each vertex x at generation m G N, we associate a sequence 
xi ... Xm where Xi G N and to simplify we write x = xi ... Xm- 

This sequence gives the complete “genealogy” of x: if y = xi ... x^ with |y| = i < m, y is the 
ancestor of x at generation i and we write y < x. 

Note that a z-ary tree is just a particular case of a GW with = 1- 

In [T], Benjamin and Lima consider the spread of an infection in a binary tree denoted T„ of height 
n. More precisely, they consider a probability vector p = (pi,..., Pfe+i) G (M!^)^+^ satisfying 

fe+i 

= i (1-1) 

i=l 

and each of the nodes of T„ is infected or not by one of the disease {1,..., fc} with the following 
rules 
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• Independently of the others, each leaf is infected according to p 

P(leaf is infected by i) = Pi, P(the leaf is not infected) = Pfc+i- 

• Nodes at generation n — 1 are infected this way: 

(Rl) if both children have the same state (infected or not), the ancestor is infected (or not) 
by it; 

(R2) if both children are infected by different diseases, the ancestor is not infected; 

(R3) if only one of the children is infected, the ancestor is infected by it. 

• This step is repeated for level n — 2 and so on. 

One of their results is the asymptotic behavior of p(n), the distribution of the state of the root of 
T„, i.e. the asymptotic behavior of 

VI < i < fc P(root is infected by i) = Pi(n), P(root is not infected ) = pk+i{n). 

They obtain the following result (0^ denotes i successive zeros) 

Theorem 1.1 Assume that Pi > P 2 > • • • > Pfe. 

If Pi = ■■■ = Pk, then p(n) converges to > 2 ^> ■ 

2. // Pi = • • • = Pi > Pi+i for some i € |1, A: — 1|, then p(n) converges to 
2 i^) • ■ •) 2 l^’ ^k-i, j where the entry 2 ^^ repeats i times. 

The aim of the present paper is to extend, when it can be, the previous results in the case of a 
z-ary tree for z > 2 and in a very specific case for a GW. Consider the family of probabilities 
define by 

( k+l 

:= < P = (pi, ■ • •, Pfc+i) e (]R+)'"+^ : ^ p, = 1 
I i=i 

and i.i.d. random vectors {0, l}^-valued such that: 

P(Xi = a) = pi,F{Xj = 1) = Pfc+i 

where (ei)(Vi are the canonical vectors of and 1 = s*- 

First we have to give the spread rules for the k diseases in a GW T„: 

• Initially each leaf x of T„ is associated to a random variable X^ and is infected as follows: 

P(a; is infected by i) = V{Xx = ef) = Pi, P(a; is not infected) = f‘{Xx = 1) = Pfc+i- 
(Consequently, each leaf is infected i.i.d. according to p.) 

• Nodes at generation n — 1 are infected this way: 

(RF) if all the children have the same state (infected or not) the ancestor is infected (or not) 
by it; 

(R2’) if two children are infected with different diseases, the ancestor is not infected; 

(R3’) if some children are infected by a single disease and the others are not infected, the 
ancestor is infected by it. 
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It can be express this way: 


for |a:| = n — 1, 


1 otherwise, 


where ® : R'" x K'^ {x, y) ^ x®y = {xiyi,X 2 y 2 , • ■ •, Xkyk)- 


• We repeat this step for level n — 2 and so on. 

As claimed, we want to determine p(n) the distribution of the state of the root and its asymptotic 

behavior, in other words the law of (or Aq) and in the case of T„: 

Vz e |l,fc],P(Ao = Ci) = Pj(n) and P(Ao = 1) = Pfc+i(n). 

In all of our results we can assume without loss of generality that Pi > P 2 > ■ ■ ■ > Pfc- 

Theorem 1.2 1. Forp S ^k such that pi = ■ ■ ■ = pk, ifp{n) converges, it does to {x,... ,x, 1— 

kx) where x is the unique fixed point in (0, i/fc] of: 

fkix) ■■= Gn{ 1 - {k- l)a;) - Gjv(l - kx) 
where Gn is the generating function of N. 

2. For p G t3^k such that pi = • • • = Pi > Pi+i > Pi+i > • • • > Pfc, then if p(n) converges, it 
does to (x,..., X, Ok-i, 1 — ix) where x is the unique fixed point of fi in (0, i/i]. 

3. If i = 1, p(n) converges to (1,0^). 

Note that the third point says that if there is only one major disease, regardless of the law of 

reproduction of N, this disease spreads a.s. to the root (asymptotically). 

In what follows, assume that N = z a.s., in other words we have a z-ary tree. 

Theorem 1.3 1. If z £ {3,4,5}, for p G tS^k such that pi = ••• = p^, p(n) converges to 

(x,..., X, 1 — kx) where x is the unique fixed point in (0, i/fc] of: 

fz,k{x) ■■= (1 - (fc - l)x)"' - (1 - kx)^. (1.2) 


2. For z G {3,4,5}, for p £ ^k such that pi = • • • = p^ > p^+i > p^+i > ••• > pfe, p(n) 
converges to (x,..., x, Ok-i, 1 — ix) where x is the unique fixed point of fz,i in (0, i/i]. 

3. If z = 6, for p G tl^k such that Pi = • • • = Pfc, x is a repelling point of fz,k- 

In section IH we study completely the case z = 6 and i = 2, where p(n) does not converge 
anymore: 

Theorem 1.4 For z = 6 and i = 2, denote 9 := {xi,Xr} where xi and Xr o,re the fixed points of 
/1 2 = fe,2 ° / 6,2 distinct of x. Then for almost all p G t^k such that Pi = P 2 > Ps > ■ ■ ■ > Pfc 

lim p(2n) = (x, x, 0 fc_ 2 ,1 — 2x) 

n—^oo 

where x £ 6. 

Remark 1.5 This result is not limited to the case z = 6 and i = 2. Indeed, denoting Xz^i = 
argmaxjQ y.-^fz,i, the only conditions a case (z,i) has to satisfy are : 

1. Xz,i is such that dxfz,i{Xz,i) < —1- 

fz,i{Xz,i) ^ Xz,i- 

3- fz,ii}h) < Xz,%- 

The paper is structured as follows: Section 2 gives the discrete dynamical system whose study leads 
to our main Theorems. The same section also gives results for the GW case. Section 3 focuses on 
2 ;-ary tree and concludes the proof of Theorem 11.31 Section 4 is devoted to the proof of Theorem 
11.41 Finally, Section 5 gives some ideas for extensions of this work. 
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2 General results for a Gallon-Watson tree 


A major part of this work consists in the study of discrete dynamical systems: given a function / 
and a value x, we study the behavior of the sequence 

In this section, we give the studied function / and some global results linked to our problem. 

To find a recursion formula, assuming that (j) is the root of its children (j)i are root nodes of 

N independent GW of height n. Then the distribution p(n+ 1) of Xq, the state of (j), is completely 
determined by the distribution of {Xi)^i, the independent states of its children with distribution 
p(n). 

Lemma 2.1 For all 1 < i < k + 1 and n > 1; 


Pi{n + 1 ) 


GAr(pfc+i(n) + Pi(n)) - GAr(pfe+i(n)). 


for i ^ k + 1 
otherwise. 


( 2 . 1 ) 


Proof : To simplify our proof denote by Sz := {1 < z < z, = 1}, the non infected sites in 
a z-sized population. According to (R1’)-(R3’): 


Pt{n + 1) = P 



2 = 2 



OO 2 — 1 




\s,\=e,<^x. 



OO 2— 1 OO 

= = '^Qz{{Pk+i{n) + Pi{n)y - 

2=2 e =0 2=2 

= GAr(pfc+i(n) + Pi(n)) - GAr(pfe+i(n)). □ 

If we define F : by : 

^ f GAr(a;fe+i + a;*) - GAr(a;fc+i), for 1 < z < fc 

1 1 - • ,a;fe+i), otherwise. 


we see that p(rz) = F'^{p). 

In fact our problem consists in the study of the fixed points of the function F. Like in [I], we first 
consider the uniform case assuming that pi = • • • = Pk- Obviously for all n > 1, Pi(n) = • • • = 
Pfc(rz) implying that we just have to study: 

Pi(rz + 1) = GAr(l - {k- l)pi(n)) - GAr(l - fcpi(rz)). 

For this purpose, define /fe : (0, i/fc] —>■ (0, i/fc] by 

fk{x) = Gn{ 1 - (k- l)a:) - GAr(l - kx). (2.2) 


We obtain the following 

Lemma 2.2 fk admits a unique fixed point in (0, i/fc]. 

Proof : Properties of generating functions ensure that fk G C°°((0, i/fe]) and on this interval: 

fl{x) = kG^il - kx) - (fc - l)G)v(l - (fc - l)a:) 

= kE[N{l - kx)^-^] -{k- 1)E[A(1 - (/c - l)a:)^"i]. 

As lim 3 ;^o+ fk{x) = E[A] > 2, fk{0) = 0, and 

fk^/k) = GN^/k) - G^(0) = ^ q,k-^ < i/fc, 

k>2 
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Figure 1: k = qs = qe = qw = i 



fk admits at least one fixed point Xk G (0, ^/k\. For x G (0, '^/k]: 

fk{x)=x O q^ ((1 - (fc - l)xY - (1 - kxY) = X 

z>2 

z-1 

qzx'^{l — {k — — kxY = x 

z>2 j=0 

z-l 

yy yy (1 — — kxy = 1. 

z>2 j=Q 


(2.3) 


X !->■ J2z>2 S^=o(1 ~ ^ -^ (1 — kxY is strictly decreasing on (0, V^] and thus bijective. 

As a result (1^ has at most one solution on (0, i/fc] implying the uniqueness oi Xk- □ 

Lemma 2.3 There exists rj > 0 such that 

'in G N, Pi(n) > a := min{Gjv(??), Pi} 

Proof : Writing Gn{Y — x) = G]v(l) — xG'm( 1) + £{x) where —>• 0, there exists r?' > 0 

ai—fO 

small enough such that if 0 < cc < ry' then |e(x)| < ^/ 2 k. Then for x < ijk =: ly and 0 < y < {k — i)x 
GAr(l - (j - l)x - 1 /) - G 7 v (1 - fx - ?/) = xG)v(l) + e((f - l)x + y) - e(fx + y) 


> xGY{l)- 


ix + y 


> G)v(l) — X > 2x — X = X, 


as G)v(l) =IE[Af] > 2 . 

If X > 77 , using the fact Gn^^ ~ Y ~ Yx — y) — GnY ~ ix — y) is decreasing in y 

K/ 1 rn 71 7^ -1 1 ~y'>~ -ix-y) 

x£[r],^/k] yG[0,(l —2x)A(fc— 2 ) 21 ] 

> inf inf GAr(l — (i — l)x — y) — GAr(l — zx — y) = inf GN{x) = GN{ji)>^- 

x^[r),'>-/k]y^[0,l-ix\ xg[7;,i/fc] 

Consequently an obvious recurrence gives Pi(n) >a = min{GA 7 ( 77 ),pij for all n > 0. 


□ 


Remark 2.1 The two previous lemmas induce that, in the uniform case, z/p(n) converges, it does 
to (xfc,..., Xfe, 1 “ kxk), which is the first point of Theorem M.Sl . 
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The following lemma ensures that the “minor” diseases can not spread to the root asymptotically 
and thus the second part or Theorem II.21 

Lemma 2.4 In the non uniform case with pi = • • • = > Pi+i, for all j > i, Pjin) —>■ 0. 

n—>oc 

Proof : Note that we just have to prove that lim„^+oo Vi+i{n) = 0. Writing Wn = ■ 


Vn > 0, Wn+i = w, 


Ez>2 9^ELd(Pi+iN+ Pfe+iW)'' ^ ^{Pk+i{n)y 


-ij=0 ' 
^z-l 


Ez >2 9^ ELo(Pi(«) + Pk+iin)y-^-^{pk+i{n)y 


< Wn 


(2.4) 


as Pi(n) > pi_|_i(n). Thus, {wn) is a positive non increasing sequence, and consequently converges. 
Denote \>y i = lim„_>+oo Wn and as wq < 1, note that I <\. 

We can find a subsequence Um such that limm^+oo Pjinm) = Oj for all j < k + 1. From lemma 
IQ we have oi > 0. Now, assume that o^+i > 0. Since f < 1, we have oi > a^+i > 0. Using (1^ 


i = liin Wn^+i 

m—>-+oo 


= lim Wn 

Tn—¥-\-oo 


^2z>2 (P*+1 (^m) H" P/c+1 (^m)) ^ (pfc+1 


?z y~) j—n (pi(rtm) T Pk+ljn-m))^ ^ ^ {Pk+lilT-m))^ 


Ez>2 9^Sj=o(«i+l +«fc+i)'' ^ -^(Qfc+i)^ 
Ez> 2 9^ Ei=d(«i + ak+iy-^-^{ak+iy 


< e 


which is a contradiction. Then Oi+i = 0 and consequently pi+i(n) —>■ 0. □ 

n—>oo 

The two previous lemmas have an important role in the following sections but we can also prove 
easily the point 3 of Theorem 1 1.2 1 According to lemma [52] lim inf „_^oo Pi(rr) > 0 and from lemma 
12.41 lim„_>oo Pjin) = 0 for all j > 1, then: 


Pi 


(n + 1) = Gat l-^pj(n) -Gjv l-^Pj(n)-pi(n) 


1=2 


1=2 


liminf pi(n + 1) = GAr(l) — Gat fl — liminf pi)?!)) = 1 — G^v fl — liminfpi(n)') . 

n—¥oo \ n—>oo / \ n—)-oo / 


Thus, liminf„_,.oo Pi(ii) is a fixed point of x i—1 1 — GAr(l — x) on (0,1]. This function being strictly 
increasing, the solution is 1. As a result hminf„_>oo Pi(ii) = liiii„_>.oo Pi(ii) = 1. 


3 Results for z-ary tree 

In this section, we investigate the special case where the GW is a z-ary tree for z > 3 (z = 2 is the 
case studied in [T]). In this case, the function F : —>• ]R^+^ is: 

F,ix) = [ (l - Sl=ioA» > forl<i<fc 

1 1 - Z)i=i Fii^)y otherwise. 

When studying the asymptotic behavior of p(n), with pi = • • • = p^ > pj+i > • • • > Pfc, it 
is enough to study the asymptotic behavior of p(n) = (pi(n), p,+i(n),... ,pfc(n)) e [0,1]'= *+U 
Indeed, we have Pi(n) = ••• = pi(n) and pfe+i(n) = 1 — X)d=iPl(^)i thus, if p(n) tends 

toward p, then p(n) tends toward p = (pi,..., pi, p 2 ,..., pfe_i+i, 1 - Ipi - Y!jZ^^ Pj). 

In the uniform case, where Pi(n) = ••• = Pfc(n), lemma [52] applies and states that there 
exists a unique fixed point ixz,k, ■ ■ ■ ,Xz,k,^ — kxz,k) of F. Furthermore, in the uniform case, 
we can restrict our study to the discrete scalar dynamical system whose dynamics is given by 
fz,kix) = Fjix, ■ ■ ■ , X, 1 — kx) for j = 1, • • • ,k and x S [0, i/fc], that is 

fz,ki^) = (1 - (fc - 1 ) 2 ;)"' - (1 - kx)"" (3.1) 
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Remark 3.1 In this section, some of the proofs use differentiation with respect to the integers z 
or k. This has to understood as a differentiation with respect to a relaxation of z or k in K.. 


The following lemma gives a lower and upper bound oi Xz k, the unique fixed point of fz k in 

(o,iA]. 

Lemma 3.1 The unique fixed point Xz,k of fz,k in (OAA] satisfies 


Hz.k ^ ^ 1; Vz ^ ^ 

where 



Proof : Recall that according to the proof of lemma [121 ^z,k{x) = fz,k{x) — x only has one 
zero Xz,k in (0, V^]- This function is positive on {0,Xz,k) and negative on {xz,k, V^] (see Figured]). 


Lower bound: We prove that, for all z > 2, and fc > 1, Az^kixz,k) > 0 which implies that 
Xz,k > Xz,k- Az,kixz,k) writes as: 


Az^k{Xz,k') — 




where = (f/zf/^’' so that Xz,k = (1 — az)lk. Obviously, Az^k{xz,k) goes to 0 as fc tends to 
infinity. Differentiating Az^k{xz,k) with respect to k: 


dAz^kix z,k) 

dk 



k-l 


(1 




(1 - az) + ^(1 


Oz) 


= ^{i-z{i-{k-\)iz,ky-^) 

< z{i-kxz,kY~y = 0 

k 

where we used the fact that z(l — kXz,kY~^ = 1- So Az^k(xz,k) is decreasing with k and its limit 
as k tends towards infinity is zero. It is thus positive for all /c > 1 and z > 2. This concludes the 
proof that Xz,k > Xz,k- 


Upper bound: Similarly to the lower bound, we prove that for all z > 2, and k > 2, Az^k{xz,k-i) 
< 0 which implies that Xz,k < Xz,k-i- As before, Az^kixz,k-i) goes to 0 as fc tends to infinity. 
Differentiating Az^kixz,k-i) with respect to k: 

dAz^k {Xz^k—Y 

dk 


> 

So Az^kixz,k-i) is increasing and its limit, as k tends towards infinity is 0. We conclude that 
Az^kixz,k-i) is negative for all z > 2 and k > 2. And in turn we obtain Xz,k < D 

In addition, to the framing of the fixed point of fz,k, the proof of lemma dUl gives that fz,k 
has a unique critical point (maximum) in (0,Yfc]. We denote by Xz,k this maximum, a direct 


_^(l_^(I_^^)z-A (1 a.) , (1 az) 


k-l 


(fc- 1)2 (fc- 1)2 


(1 - az) 
(fc- 1)2 

(1 - aY 

(fc- 1)2 




(1 - z(l - kxz,k)'' ) = 0. 
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computation of the critical point of fz^k in (0, V^] yields the value (ld.2l) . Furthermore, another 
direct computation of the zeros of d^fz,k gives the unique inflection point x* f. of fz,k in (0, i/fe]: 




- {k - 1)^-1 

- (fc- 1)^’ 


_ fc ^-2 - ( fc - 1)^-2 
k^ - {k - 1)^ 


(3.2) 


Note that Xz,k < x*k- 

Now, we show that having Xx,k as an attracting fixed point of fz,k guaranties that the asymptotic 
behavior of the diseases spread is as stated in our main result. 

Proposition 3.2 For the uniform case, if Xz,k is such that \dxfz,k{xz,k)\ < then p(n) converges 
to {Xz,k-i ■ • ■ 5 ^z,k^ 1 kXx^kf 

Proof : Since our result does not depend on z and k, we will drop the indexation of /, x and x 
by those integers. 

Note that \ix>x, necessarily f'{x) S [0,1) and /((O, i/fe]) C (0,/(x)] C (0,x]. Since f{{0,x]) = 
(0, x\ and f{x) > a: on this interval, /"(x) tends toward x. This result is still true if x G [x, x] with 
a similar reasoning. 

In the rest of the proof we assume that x < x and it is stuctured as follows. First, we study the 
variations of /^ = / o / and then show that only has one fixed point. Finally we conclude on 
the convergence of the discrete scalar dynamical system. 

Variations and concavity of 

First note that as x < x, /(x) > x and since / is increasing on [0, x] there exists a unique xi G (0, x) 
such that f{xi) = x. Similarly, / is decreasing in [x, V^]) if /(V'^) ^ ^ there exists Xr G (x, V^] 
such that /(xr) = x and none otherwise. 

Thus, as (/^)'(x) = /'(/(x))/'(x) and /([O, i/fc]) = [0, /(x)] C [0, Vfc], the only critical points of 
are xe, x and possibly Xr- Since /^(x^) = /(x), it is a maximum of (in [0, V^])) liks Xr in case of 
existence, and thus x is a local minimum. To conclude, is increasing on [0, Xi] U [x, min(i/fc, Xr)] 
and decreasing elsewhere. In particular, is increasing on [x,x] C [x, min(i/fe, x^)] (we assume 
that if Xr does not exist in [0, i/fc], we take Xr = +oo). 

Note that this study of the variations of is not restricted to the attracting fixed point case, as 
long as we have x < x. 

A study of the second derivative of p gives that it is negative on [0, xi] and cancels only once 
on [x,Xr.]. So p is concave on [0,x^] and have exactly one inflection point on [x,Xi.] where it is 
convex then concave. Note that this study of the inflection points of p is inconclusive on [x^jX^] 
(with Xi the first pre-image of x* by /) and [x*, ^/k\. 

X is the unique fixed point of p in (0, i/fc]: 

We reason by contradiction with the maximum number of inflection points of p in [0,Xf] and 
[x,Xr]. 

First note that as x is such that /'(x) G (—1,1), (PYP) G (0,1) and so there exists ry > 0 such 
that P{x) > X if X G {x — f],x) and P{x) < x if x G (x, x + ry). 

Now, assume xi = supjx G (0,x),/^(x) = x} exists (necessarily xi x), thus p{x) > x on 
(xi,x) and (/^)'(xi) > 1. We face two cases: either xi > x, either xi < x (as xi = x gives directly 
a contradiction since {p)'{x) = 0). 

Assume xi > x: as (/^)'(xi) > 1, {P)'{x) = 0 and {P)'{x) < 1, P has at least one inflection 
point in (x,x). Since P{xi) = xi, then X 2 = /(xi) = inf{x G {x,Xr), P{x) = x} is also a fixed 
point of p. Similarly as before, p has at least one inflection point in (x,Xr). So p has a least 2 
inflection points on [x, x^.] which raises a contradiction. 

Since p has no fixed point in [x,x), it has none in [x^,x] because p is decreasing on this 
interval. Moreover, p is concave on [0,X£], P{0) = 0 and P{xi) > xi so that p{x) > x on {Q,X(\. 

We thus conclude that xi does not exist. If there exists a fixed point X 2 of p in (x, i/fc] then 
/(X 2 ) G (0,x) is also a fixed point of p which is not possible. So p doesn’t have a fixed point in 
(x, i/fc] either. 




Remark 3.2 Note that a corollary to the uniqueness of the fixed point x is that f‘^{x) > x. 

X is asymptotically stable for f in (0, ^/k\'. 

We assume that Xr exists, otherwise the reasoning is the same with Xr replaced by i/fc and one less 
step: we would not have to consider the interval [xr^ V^]- 

From the variations of /^, we have that f'^{[x,x]) = [f^{x),x] C So if a: £ [x^x] then 

/^"(cc) is an increasing bounded above sequence, and so it tends towards x, the unique fixed point 
of p in [x,x\. As X is the fixed point of / and / is continuous, / o p'^{x) = also tends 

towards x. 

Since f^(xr) = f{x) < Xr, we have that f'^{[x,Xr]) = [x,f{x)] C [x,Xr]- So if x £ [x,Xr], then 
/^"(x) is a decreasing bounded below sequence. As for the [x,x] case, we conclude that /"(x) 
tends toward x. 

As p{x) > X and P{x() = f{x) < Xr, we have that p{[xi,x]) = [f^{x), fix)] C [x,Xr]. And 
again, we conclude that if x £ [x£,x], then /"(x) tends towards x. 

We have that /^((O, x^]) = (0, /(x)] and if x £ (0, x^], then /^(x) > x so /^"(x) is an increasing 
bounded above sequence so it converges to x. Finally, /^([xr, V*^]) = [/^(V^); /(^)] C (0, /(x)] and 
the previous case gives the convergence to x. 

So we have that for all x £ (0, ^/k], lim„^oo f^{x) = x, which concludes the proof. □ 

Figure [5] sums-up what we gather of the variations of fz,k when x^^k is linearly attracting. 


Figure 2: Table of variations of fz,k when Xz,k < Xz,k 



Proposition 3.3 If Xz,i is a linear attractor for the uniform case with i diseases, then x = 
ixz,i ,..., Xz,i, 0fc_i, 1 — iXz,i) is a linear attractor for the non uniform case with Pi = • • • = Pi > 
Pi+i > • • ■ > Pfe- 


Proof : X is a linear attractor if all the eigenvalues of the matrix A := §|-(y) are in (—1,1), 
with F = (Fi, • ■ ■, Fk-i+i) = {Fi,Fi+i, ..., Fk), a truncated version of F and y = ixz,i, Ok-i)- F 
is defined by 


(l - (z - l)xi - Xj) - if ^ = 1 

(^1 - ixi - Xj'j otherwise. 
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Note that ^f(x) equals to 


-z{i - 1) (l - (i - l)a;i - ^o) + (l - ixi - ^j) 

-iz (l - ixi - ^j) +iz{l- ixi - Xj'j 

z (l - ixi - Xj'j 

-z(l- iXi - Ej=2j^ ^j) +z{l-iXi- Xj'^ 


a m = i = 1 
it m 1, £ = 1 
it m = i > 1 
itl^m^i^l 


implying 


dPrn 

dxe 


(y) 


-z(z - 1)(1 - (i - l)a;^,i)"= '^ + zi{l - ^ 

0 

z(l - ix^^iY-'^ 

0 


if TO = £ = 1 
if TO 1, £ = 1 
if TO = £ > 1 
ifl^^TOT^f^l 


Thus A is upper triangular and its spectrum is {—z{i — 1)(1 — (* — l)xz,i)^ ^ ~ ^^z{l — 

iXzY^~^}- First note that as Xz,i is an attracting fixed point for the uniform case then ff^(y) € 

(- 1 , 1 )- 

It remains to prove that (for j = 2,... ,k — i + 1), 


0 JY ■ 

0 < -^ixz,i,0k-i) = z{l - iXz,iY~^ < 

and note that positivity is obvious. 

From lemma EH Xz,i > Xz,i for all z > 2 and i > 1, and 


OF _ dF 

~(.Xz^itOk—i') ^ (^Xz.itOk—i') — 1 , 


dxi 


which concludes the proof. 


□ 


Lemma 3.4 For the uniform case, and z > 2, there exists Kz > 0 such that, for all k > Kz, the 
fixed point {xz,k, ■ ■ ■, Xz,k, 1 — kxz,k) is attracting if z < 5, repelling otherwise. 

Proof : It is enough to study the asymptotic behavior of dxfz,k{xz,k), the non zero eigenvalue 
of the linearized dynamical system. 

From lemma o and the continuity of dxfz.k'- 

lim dxfz,kixz,k) = lim dxfz,k{xz,k) 

fe—foo k—>oc 


Furthermore, the expression oi Xz^k gives us an equivalency as k tends to infinity 


dxfz,k{S^z,k^ 


-z{k - 1)(1 - (fc - l)xz,kY~^ + - kxz,kY~^ 

-z{k - 1)(1 - (fc - l)xz,kY~^ + k 

k{l - z(I - {k- l)xz,kY~^) + z:{t - {k- l)xz,kY~^ 


k—^oo 


—y 

k—¥oo 


k 


^1 - z{k - kxz,kY ^ 



^z.k 


1 — kxz 


I - (z- l)z^ 



+ 1 


And this limit is in (—1,1) for z < 5 and is less than — 1 for z > 6, which concludes the proof. □ 
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Remark 3.3 If we were able to prove that the eigenvalue dxfz,k{xz,k) is decreasing with respect 
to k, then its asymptotic value would directly give us that for z < ^, the fixed point is a linear 
attractor for all k > 1. Even though this monotonicity seems true numerically, we were not able 
to prove it, and were reduced to cumbersome computations for z € {5, 6 }. 


The following result focuses on the fixed points of 3-ary, 4-ary and 5-ary trees: 

Lemma 3.5 For the uniform case, if z € {3,4,5}, and k > 2 the fixed point Xz,k is linearly 
attracting, that is dxfz,k{xz,k) S (—1,1). 


Proof : For the upper bound on the derivative, we have, for all z > 2 

dxfz,k{xz,k) = -z{k-l)z{l - {k - l)xz,ky~^ + kz{l - kxz,ky~^ 

< -z{k - 1)(1 - kxz,ky~^ + kz{l - kxz,ky~^ 

< z(l kXz,k) 

< z{l - kxz,ky~^ = 


the last inequality comes from lemma [XT] 

For the lower bound on the derivative, we treat differently the cases z G {3,4} and z = 5. 
Lower bound, case z G {3,4}: noting that the unique inflection point of fz,k is such that 

^z,k ^ ^z,k 

dxfz,k{x% y = min dxfz,k{x)<dxfz,k{xz,k) 

’ a:e[ 0 .Vfc] 

Direct computations for z = 3 and 4 give 


dxh,k{xlk) = 


3(fc — l)k 
3fc2 - 3A; + 1 


^ 1 ) fA,kix^y) 


4(fc — l)k 

(1 - 2fc)2 


> -1, Vfc > 2, 


which concludes this case. 

Lower bound, case z = 5: here dxf 5 ,k{x^y < —1 for /c > 3 so we cannot use the same 
argument. Instead, we prove that x^^k-i < x^ j. for k large enough, which give dxf 5 ,k{x 5 ,k-i) as a 
lower bound for dxf 5 ,kix 5 ,k)- We write 


and 




fik 

1 + kfik 


Pk 


1 


5 

k-V 



- 1 > 0 


1 

4 


> 0 


Comparing x\ ^ and x^^k-i 


- Pk{k-l-dk)-5 

X,^k-X,,k-l- 

The denominator of which is positive. If fc G C is a zero of the numerator then, straightforward 
computations give 

0 = k^{k - 1 - 5kf -{k- 1)2(<5 + (fc - 1 - 5k)f G M4[^] 


A numerical computation of the roots of this fourth order polynomial gives 2 complex conjugated 
roots and 2 real roots. The 2 real roots are approximatively ki = 0.3079371 and k 2 = 3.3623924. 
So the numerator has constant sign for k > k 2 and in particular for fc > 4. Noting that 2:5 4 > x^^s 
we conclude that x^ j^> x^^k-i for all fc > 4 and the following lower bound for the eigenvalue 

dxf5,k{x5,k) > dxfb,k{Xb,k-l)^ V/c > 4 
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A straightforward computation gives that dxf 5 ,k+ 27 ix 4 ^k-i+ 27 ) > —1 is equivalent to 


(6 - 4 • + (634 + 6\/5 - 432 • + (25126 - 17496 • 5^/'‘486\/5 - 4 • 

+ (442634 - 314925 • + 13122\/5 - 216 • b^/^)k 

+2924642 - 2125764 • + 118098\/5 - 2916 • 5^^^ > 0, 

which is true for all A: > 0. So dxfh,k{S:b,k-i) > — 1 for all k > 27. 

For A: G |1, 26], a direct numerical computation of the eigenvalues gives that they all belong to 
(-1,0). □ 

Remark 3.4 It is actually possible to show that dxfz,k(xz,k) < 0 for z G {3,4,5}. To do this, 

one can compute dxfz,k{xz,k) cind show that it is decreasing with respect to k > 3 for z = 3 

and decreasing for k > 2 for z G {4,5}. This leads to dxfz,k{xz,k) < 0 for these {z,k) and as 
Xz,k < Xz,k, Xz,k > Xz,k- So that fz,k is decreasing at Xz,k- 

Lemma 3.6 For the uniform case, z = 6 and k > 2 the fixed point xe^k is repelling. 

Proof : We prove that dxfe.ki^e.k) < -1, using the fact that dxf6,kixe,k) < 7aa7K{dxfe,k{xz,k), 
^xfe,k {Xzjk—l^f 

As for the case z = 5, dxfe.k+iixe.k+i) < —1 is equivalent to 

(7 - 5 • 6^/^)A:® + (125 - 75 • 6^/^ - 10 • + (900 - 450 • 6^/^ - 120 • - 10 • 

+ (3265 - 1350 • 6^/® - 540 • 6^/® - 90 • 6^/® - 5 • 

+ (5960 - 2025 • 6^/^ - 1080 • 6^/® - 270 • 6^/® - 30 • 6^/®)A: 

+ (4373 - 1215 • 6^/® - 810 • 6^/® - 270 • 6^/® - 45 • 6*^^) < 0, 

which is true for A; > 0. So dxfe,k{xe,k) < —1 for all A: > 4. 

A similar computation for dxfe.k+iixe.k+s) leads to dxfe.kixe.k-i) < —1 for all A; > 4. So 
dxf6,k{xe,k) < -1 for A: > 4. 

A numerical computation of dxfe.kixe^k) for k = 2 and 3 yields values less than —1, which 
concludes the proof. □ 


Remark 3.5 Note that as dxfix*f.) > —1 for z G {3,4}, the fonction fz,k is contracting on 
[xz,k, V'=]- This gives the main ingredient to an easier way to prove the convergence of the uniform 
case than in Provosition 1 3 . ill However, the given proof of Proposition\K^is more general as it only 
requires the fixed point to be linearly attracting. 

Lemma 13.51 in addition to Proposition 13.21 gives the convergence in the uniform case for z G 
{3,4, 5}. The following Proposition extends this result to the non uniform case. 

Proposition 3.7 Assume that, for the uniform case of the z-ary tree, p(n) converges to {xz,k, ■ ■ ■, 
Xz,kA ~ kxz,k)- Then, in the non uniform case with pi = ••• = Pi > P*+i > • • • > Pfc, P{n) 
converges to {xz,i,..., Xz,i,0k-i, 1 - ixz,t)- 

Proof : Assume that in the uniform case p(n) converges. Obviously if pi = • • • = Pi and 
Pj = 0 otherwise, p(n) converges to {xz,i, ■. ■ ,Xz,i,0k-i, 1 — ixz,i), where Xz,i is the fixed point of 
fz,i- We want to extend this result to the case pi = Pi > Pi+i > • • • > Pk. 

We introduce the function F = {Fi,..., Fk-i+i) = {Fi, Fi+i,..., Fk) which is a truncated 
version of F 

((i - (* - i)xi - ew "' - (i - 
\ (l - xj'^ ~y~ Xj^ 
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if ^ = 1 

if ^ = 2, • • • ,k — i + 1, 


defined in the set i := {{xi,... ,Xk-i+i) S > X 2 > X 3 > ■ ■ ■ > Xk-i+i > 0,ixi + 

< 1}. So, equivalently to our convergence result, we show that, for x G ^k,i, i^”(x) 
converges to Xi = {xz,i, 0 k-i)- 

Let e > 0 and p = (pi,... ,pfc_i+i) G ^k,i- From lemma 1^31 we have that, for all n G N, 
Pi(n) G [a, V*]- Now, while noting that F'^{x,0k-i) = (/"j(a;), Ofc_i), we define En = {x G 
[a, ^/i],fz,i{x) G B{xz, 1 ,^/ 2 )}. Clearly, from the convergence of the uniform case 

[a, V'] C U„>o£ln- 

As the inverse image of an open set by a continuous fonction, En is also an open set for all n G N. 
Since Un>oEn is a sequence of open sets covering the compact [a, i/i], there exists N such that 

[a, V*] C 

implying that Wx G [a, V^]) F^ix,0k-i) G B{xz, x {0}'"“* C B(xi,e/ 2 ). 

On the closed set := [a, i/i] x fl ^k,i, is uniformly continuous and thus there exists 
i5 > 0 such that 

y{x,y), (x',y') G Ux,y) - {x',y')\\ <6^ \\F^{x,y) - F^{x',y')\\ < ^2 

According to lemma [231 Pj{n) = F"(p) —^ 0,Vj > i, and consequently there exists A^i G N such 
that if n > A^i, ||(pi+i(n),... ,pfc(n))|| < S implying 

||F^(pi(n),p,+i(n), ... ,Pfe(n)) - F^(pi(n),0fc_i)|| < e/2. 

Thus, recalling that according to lemma [231 Pi(u) G [a, V*] 

||(pi(n +Af),p*+i(n +Af),...,pfc(n +AT)) -x|| = ||F^(pi(n), p*+i(n),..., pfc(n)) - x|| 

< l|F^(PiH,P*+iH, ■ • ■,Pkin)) - F^(pi(n),0fc_i)|| + ||F^(pi(n), 0fc_i) - x|| < e/2 + e/ 2 = e. 


□ 


4 A Specific case of Repelling point and attracting orbit 

In all this section we study the case of a 6-ary tree with 2 dominant diseases. According to lemma 
na the convergence to the fixed point is no longer true and in the present section we prove the 
existence and uniqueness of an attracting orbit of prime period 2. 

For the sake of clarity, we write x instead of xq ^2 and likewise for /, F, x*, x, x. Moreover, an easy 
fact is the following 

Remark 4.1 If y is a fixed of but not of f, then f{y) is a fixed point of f^ distinct of y and 

{fYiy) = (fnfiy))- (4.1) 


Lemma 4.1 All the fixed points of f'^ are in (x,Xr) where Xr '■= sup{x G [0, Y 2 ],/(x) = x} and 
recall that x = argmaxjg iy 2 ]/(x). 

Proof : First note that x < x and: 


fix) > 0 X < X. (4.2) 

As X = argmax2,g[o,i/2] fix), /([0, 1 / 2 ]) = [0,/(x)] implying that x G [0,/(x)] and the existence of 
XI G [0,x) such that fixf = x and fixe) = fix) = max,j,g[oy/2] fix). 

Moreover note that 

f"ix) > 0 X > X*, 
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Figure 3: z = 6,i = 2, f and p 



Figure 4: Table of variations of 



X* is the only inflexion point of / implying that /' is decreasing on [0, a;] as a;* > a;. 

First we show that there is no fixed point on (0, a:) and we split this in two cases: (0, X(\ and [x(,,x). 
On (0,if]: 

Using ■ra, on this interval: 

sgn{f{x) - fix))' = sgn{f'{x)if'{f{x)) - 1)) = sgn(/'(/(a:)) - 1). (4.3) 

/' being decreasing on [0,x], /'(/(a:)) — 1 also on [0,xe), and as /'(/(O)) — 1 = /'(O) — 1 = 5 and 
f'ifixf)) — 1 = fix) — 1 = —1, there exists /3 S (0, xf) such that f — f is non decreasing on (0, /3) 
and decreasing otherwise. Since /^(O) — /(O) = 0 and fixf — fixf = fix) — x > 0, fix) > fix) 
on [0,5^]. Thus, as fix) > x on (0,a;^), f do not admit a fixed point on this interval. 

On ixi, x): 

fix) > 0 and as fix) e (x,/(x)), fifix)) < 0 : 

(fix) - x)' = fix) fifix)) - 1 < 0. 

Thus fix) — x is decreasing on fi, x) and, as algebraic computations give fix) — x > 0, fix) = x 
admits no solution on (4^,4;]. 

In a second time we make a similar reasoning on [xr , V^] ■ The existence of Xr is clear noting that 
/( 1 / 2 ) = 1 / 2 ® < X. 
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On [Xr, V2): 

f'{x) < 0 and /([x^, 1/2]) = [/(V2),i]- For all y € (/(V2),x), f'{y) > 0, thus: 

Vx e [xr,y2),{f)'{x) < 0. 

As a result p{x) — x is decreasing on [x^, V2) and as < 1/2, P has a unique fixed point on 

this interval if and only if P{xr) > Xr- 

Assume that /(Xr) > Xr and denote by j/ G [Xr, 1/2) a fixed point of /^. According to remark HTTl 
f{y) is also a fixed point of p (not of /) and 

{PPPy)) = {P)\y) < 0. (4.4) 

But as p is increasing on (0,x^) U (x,Xr) and not increasing elsewhere, (14.411 implies that f{y) G 
(x(,x), which is a contradiction as we have proved that there is no fixed point on this interval. 
Consequently all the fixed points of p are on (x,Xr)- □ 

Figure 5: z = 6,i = 2 zoom on the fixed points of p 



Lemma 4.2 There exists exactly two fixed points of p distinct of x. They are both attracting. 
Proof : Noting that fix) < —1, x is a repelling fixed point of p such that 

ipyp) = {fpp > 1 - 

So there exists £ > 0 such that p{x — e) < x — e and x + £ < p{x + s). As P{x) > x and 
P{xr) < Xr, there are at least two additional fixed points on (x,Xr). Moreover the fact that p is 
increasing on this interval ensures the existence of two attracting fixed point from either side of x. 
We have already seen that the fixed points of p are in [x, Xr] and moreover a study of the sign of 
{Py shows that p admits exactly one inflection point in [x,Xr]. As a result p has three fixed 

points in (0,1/2]. □ 


Proposition 4.3 In the uniform case, if X is the Lehesgue measure on K; 

A Qpi G [0,1/2], p(n) {x^Xrlj^ = 1 / 2 . 

Proof : According to the lemmas [4.11 and [4.21 p has exactly three fixed points on [x,Xr], 
denote xi and Xr respectively the one in (x,x) and the one in (x,Xr). 

On (x,Xr): 
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is increasing, f'^{{x,Xr)) = {x,Xr) and P{x) > x then /^(x„) = Xn +2 > Xn- The sequence 
{x 2 n) is increasing and bounded above by Xr, thus it tends to a fixed point of in {x,Xr], Xr- 
On [Xr, Xr]'- 

f'^{[xr,Xr]) = [xr,f{x)] C [xr,Xr)- Moreover on this interval, f^{x) < x implying that X 2 n is 
decreasing and bounded below by Xr, thus convergent to Xr- 
On {x, x): 

the reasoning is the same as the one on {x, Xr), and for all x on this interval X 2 n xi- 
On {xi,x): 

there exists an unique xi such that f{xi) = x. Note that on this interval p is decreasing implying 
that 

• if a: G [x^, xi), p{x) G (S, f{x)] C (x, Xr) and X2n Xr', 

• if a; G (a;i,a:], P{x) G [p{x),x) C {x,x) and then X 2 n —>■ Xi- 
On [xr^ 

there exists X 2 such that f(x 2 ) = x. On this interval p is decreasing and 

P{[Xr,X 2 )) = {x,f{x)] and p{{x 2 ,p]) = [p{y 2 ),x) 

If a; G [xr, X 2 ), P(x) G (a;, f{x)] and we are in one of the previous cases. Let us postpone the other 
case {x G {x 2 , V^])- 
On [0,Xf]: 

P is increasing implying the existence of 7 such that pp) = xe- If a; G [ 7 , 3 ;^], P{x) G [xi,f{x)] 
and we are in one of the previous cases. 

We can easily show by contradiction that for all x G [ 0 , 7 ] there exists n such that P'^{x) G 
[x£,f{x)]. If it is not the case P^{x) is an increasing sequence bounded above, and as a conse¬ 
quence it converges to a fixed point of p in [0, xi]. According to lemma im there is no such point. 
There is however countably many points that converges to x (because /^([0,7]) C [0,/(a;)]). 

To conclude if a; G ( 3 ) 2 , V^]; P{x) G [p{^/ 2 ),x) and we are in one of the previous cases. □ 


To obtain results for the non uniform case, we need the following: 

Lemma 4.4 {(a:^,a:^, 0 fc_ 2 ,1 — ^xi), {xr,Xr,0k-2, 1 — 2xr)} is an attracting orbit for F in 

Proof : Note that it is the same to prove that {{x£,0k-2), {xrPk- 2 )} is an attracting orbit 
for F = (Ai,..., Fk-i) where 


Feix) = 


(1 - xi - J2j=2 ^j) ~ ^j) 

(1 - 2xi - Xj^ ~ 


if£ = l 
otherwise. 


Let Bij(x) = 


dxi 


“ Em=l dx„ 


{F{x))^^{x) and note that: 


BiJp£y0k—2) — Bijprj^k — 2 ') — ^ 




, vk— 2 ) otnerwise. 


^ixi,0k-2)^iXr,0k-2) 


< I accord- 


Ok- 2 ) is upper triangular and 0 fc_ 2 )| = 

ing to lemma It remains to prove that the other eigenvalues of B are (strictly) bounded above 
by I, the positivity being obvious. 

With the same notations as the one of the proof of lemma IXTl recall that if x < i, §|^(x, 0 ^- 2 ) < 1 
and greater than I otherwise, and one can easily check that P{x) > x. Then, as x < x, obviously 
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X <xt and consequently ^^{xi,Qk- 2 ) < 1 and as i < Xr, 0 /c_ 2 )f|^(a;r, 0 fc_ 2 ) < 1 . □ 

Now we are able to prove Theorem 11.41 The previous lemma ensures the existence of e > 0 
such that if X S B{{xr,Qk- 2 ), s) (respectively B{{xg, Ok- 2 ), £)), then lim„_j.+oo F^"(x) = (x^, 0 ^- 2 ) 
(respectively (x^, 0 fc_ 2 )). 

According to the stable manifold Theorem, the stable manifold of F, lT®(x, 0 ^- 2 ) is a one di¬ 
mensional smooth manifold in a neighborhood of {x, 0 k- 2 )- So there exists e' such that W^, = 
lT®(x, 0 fe_ 2 ) n i?((x, 0 fc_ 2 ),e') > 0 is negligible with respect to the Lebesgue measure. In other 
words, for almost every x £ B{{x, Ok- 2 ), s'), there exists m > 0 such that F^{x) ^ B{{x, Ok- 2 ), s')- 
Let En = {xi G [x -I- £', Xr], .F^"(xi, 0 fe_ 2 ) £ i3((xr., 0 fc_ 2 ), ®/ 2 )}, and according to Proposition 14.31 
we have already seen that [x-be', Xr] C U„>oA„. With a similar reasoning as the one of Proposition 
[pi there exists N > 0 such that [x + e',Xr] C U^^QEn implying that for all xi £ [x -b s',Xr\, 

F2A^((xi,0fc_2)) £ B{{Xr, 0 k- 2 ),^l 2 ). _ 

Using the uniform continuity of on the compact set ■.= [x + s',Xr\ x n i^k, 2 , there 
exists 12 > 0 such that 


Vx,y £ ||x - 2/11 < ^ ^ ||F^^(x) - F^^( 2 /)|| < -/ 2 . 


Thus for X = (xi,. .., Xk-i) £ [x -b s', x^] x B( 0 k- 2 , v)\ 

||F^^(xi,... ,x/c_i) - (Xr,0fc_2)|| < ||F^^(xi,... ,X/C_i) - .F^^(xi,0fc_2)|| 

< e /2 -b e /2 = e. 

Consequently F^^(x) £ B{{xr,0k-2),s) implying that F^"(x) —>■ Xr- 
Note that we can realize exactly the same reasoning on [x, x — e']. 

Let xi £ [a, V2]\[x, Xr] (see Lemma for the definition of a) and according to Proposition 14.31 
there exists n such that F^"(xi,0fc_i) £ [x,Xr] x i?(0fc_2,i^). Using again the uniform continuity 
of F^", there exists v' such that 


X - : 


<v' ^ ||F2”(x) - F2”(y)|| <5:= min v. 


F/"(xi,0fc_2)-X Xr-F/"(xi,0fc_2) 


implying that if |l(x 2 ,... ,Xfc_i)|| < v' then 

X + F/"(xi,0fe_2) ^ ^ ^ Xr +F/"(xi,0fc_2) 

2 S xi [xi,... ,Xk-i) < 2 

0< F|”(xi,... ,Xfe_i) < iz,Vj > 1 
and thus F^"(xi,..., Xk-i) £ [x, Xr] x B{0k-2,v). 

Note that we write p = (pi,..., pfe_i) for (pi, Pa, P 4 , ■ • ■, Pfe) and p(n) = F"(p). 

According to lemma lOl Pj(n) —>■ 0 for all j >2 and according to lemma H31 Pi(n) £ [a, 1 / 2 ] for 
all n > 1. There exists Ni such that 'in > Ni, p(n) £ ( 0 , 1 / 2 ] x B{0k-2,v'). Then there exists 
n > 0 such that 

p(iVi + 2n) £ [x, Xr] X B{0k-2, v). 

If p(A^i + 2n) £ [x,x — e'] U [x -b s',Xr\ x B{0k-2,Js), we have the convergence to (x^,0fc-2) or 
(xr,0fe_2). Otherwise p(Aii +2n) £ B{x,e') x B(0k-2,v) and 

• either lim„^+oo p(2n) = (x, 0 fc_ 2 ), 

• or there exists m > n such that p(Aii + 2m) £ [x, x — s'] U [x + s', Xr] x F(0fe_2, v) and we 
have convergence to (x^, 0 fe_ 2 ) or {xr, 0 k- 2 )- 
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If we are in the first case, it means that p G = U„>o-F and as a countable union of 

negligible sets with respect to the Lebesgue measure, Wg is also negligible: 

A(pG.^fe, 2 , lim p(2n) ^ {(s^Ofc- 2 ), (ir,0fc-2)} ) = 0 (4.5) 

\ n —>^+00 / 

Remark 4.2 This result is not limited to the case z = 6 and i = 2. Indeed, according to our proof, 
the only conditions a case {z,i) has to satisfy are : 

1 . Xz,i is such that dxfz,i{Xz,i) < —1- 
fz,i{flz,i) ^ d^zp- 

3- fz,i{^l^) < Xz,i- 

The last two conditions appear in the proof of lemmas an d \T^ While condition 3 seems to be 
always satisfied, condition 2 is not. For instance, it is not satisfied for (z, i) = (12, 2) and numerical 
simulations show the existence of an orbit of period 4 (see Figure\^ for this case. However, the 
aim of the paper is not the study of all non convergent cases. 

Figure 6: z = 12,i = 2, / (red), (yellow), f^ (blue) 



5 Open questions and variant case 


In [I], the authors study a binary tree with the following changes 
• k = 2, we have only two diseases; 


• (R3) is replaced by (R4): if only one of the leaves is infected, the node is infected by it with 
probability a and not infected with probability 1 — a. 


The authors obtain the following 


Theorem 5.1 For all p G 13^2, Pin) converges and 


1 . Ifa> 1/2 


lim p(n) 

n—>-oo 


( 1 , 0 , 0 ) 

( 2ol-1 


20-1 1 \ 
4a:— 1 ’ 4a— 1 j 


if Pi > P2, 
*/Pl = P2- 


2. If a = 1/2 


lim p(n) 

n—^oo 


(Pl - P2,0, 1 - Pi +P 2 ) *fPl>P2, 
(0,0,1) i/Pi=P2- 
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3. Ifa< 1/2 lim,i^oo p(n) = (0,0,1). 


In the same spirit, we replace (R3’) by (R4’): 

where Ui is the probability of infection i and ^i{z) := {S^ = {1 ,..., = ■ ■ ■ = = Cj} 

the event the I firsts children are not sick and the others are infected by disease i. 

Recall (Rl’)- if all the children have the same state (infected or not) the ancestor is infected (or 
not) by it (which does not seem natural). 

Figure 7: (R4’) for = 3 a.s. 


= {1-a.r 



We obtain the equivalent of lemma 12.11 
Lemma 5.1 For all n > !.• 


Pi(n+1) = 


G'Ar(pfc+i(n) + Pi(n)) - GN{v>k+i{n) + p*(n)(l - ai)) + GnU^ - ai)pi{n)), if i < k 


1 — Pj(’2- + 1) otherwise. 

Proof : For i ^ k + 1, the recursion formula obtained here is: 

00 Z — 1 


p.,(n + 1) = ^ ^ P Xj = , 


2=2 ^=0 

00 2— 1 






)Xj = a \£/f{z) 


2=2 ^=0 
00 f-S —1 


vi=i 




2 = 2 


^(n)(l - (1 - Oi)'' ^) + Pi(n)^ 


.e=i 


= ((Pfe+i(’^) + Vi{n)y - (pfc+i(n) + Pi(n)(l - ai)Y + (1 - a*)"'pi(n)^) 


= G'Ar(pfc+i(n) + Pi(n)) - G'Ar(pfe+i(n) + Pj(n)(l - o^)) + G'iv((l - ai)p^{n)). □ 
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In all the following we assume that 


Vie ll,k},ai=a>E[N]-\ 


(5.1) 


Lemma 5.2 There exists rj > 0 such that 

Vn e N, pi(n) > min{GAr(7?),pi} 

Proof : According to (j5.1ll . there exists d > 0 such that a = + 5. As Gm{^ — x) = 

Gtv)!) — a:G(Y(l) + s{x) where ^G)/x 0 when x goes to 0 , there exists ?? > 0 small enough such 

that a 0 < X < rj, |e(a:)| < Then for 0 < a; < ^/k and 0 < y < {k — i)x: 

Gj\r (1 - (i - l)a; - y) - Gn (1 - (i - 1 + a)x - y) > aa;G(v(l) - 6 xG'pf{l) > x, 

and as GAr((l — a)x) > 0 : 

G 7 v (1 — {i — l)a; — y) — Gn{^ — (i — 1 + a)x — y) + GAr((l — a)x) > x. 

The rest of the proof is exactly the same as the one of lemma 12.31 □ 

Lemma 5.3 For all j > i, PiM —>■ 0. 

n—^oo 

Proof : Writing for all 1 < J < fc: 

Z — 1 

Pjin + l) = a'^q^pj{n)'^{pj{n) + pk+i{n)y~^~^{pk+i{n) + (1 - a)pjY + (1 - a)''p*+i(n)^ 

z>2 £=0 


it is not difficult to see that the sequence w'„ = is positive and decreasing. The rest of the 

proof is the same as the one of Lemma 12.41 □ 

Proposition 5.4 //pi > P 2 > ■ ■ ■ > Pfe then lim„^+oo Pi(?^) = 1- 

Proof : We have already shown that p, (n) —^ 0 for all j > 1 and according to lemma 15.21 

n—^'OO 

liminf„^oo Pi(?T-) > 0. Consequently 


k 


k 


Pi(n + 1) = Gat l-^Pj(n) -Gat l-^Pj(n)-api(u) +Gjv((l-a)pi(n)) 


liminf pi(n + 1 ) 

n—>-oo 


i=2 


i=2 


Ga/(1) — Gn — aliminf Pi(n)') + Gn f(l — a) liminf Pi(u)') 

\ n—^oo J \ n—voo / 

Gn (1 — aliminf pi(n)) + Gn ((1 — a) liminf pi(n)) . 

\ n—¥oo / \ n—¥oo / 


1 - 


Thus, liminf„^oo Pi('ra) is a fixed point of x 1 —>■ 1 — Gat)! — ax) + Gat))! — a)x) on (0,1]. This 
function being increasing, the solution is 1. As a result liminf„_).oo Pi(* 2 -) = liiii„_>.oo Pi (''2-) = 1- 

□ 


Here we give some open questions that may be interesting to study. 

1. What happens if we slightly change (Rl’) and (R2’): although all children are sick, we need 
at least a disease to be expressed. In a more probabilistic way, for all 0 < f < z — 1 



— Ci 


^Hz)\ = 


W — 1 
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Assuming that ai < E[A] 


(5.2) 


lim Pi{n) = 0. 

n—v+oo 

Indeed, we can bound above lim„_i.+oo pi(n) by the probability that in a GW with repro¬ 
duction law N and percolation probability ai, there is an infinite connected component 
containing (j). It is well known that ac = E[7V]“^, is the critical probability for the existence 
of such a connected component which gives 

2. An interesting subject is the link between the law of N and the existence of a fixed point. It 
seems obvious that if N is heavy tailed, there is no fixed point but what happens if P(A^ > 6) 
is very low? Is there a critical value for this probability implying the existence or not of a 
fixed point? 

3. Numerical simulations suggest the existence of a unique attracting orbit for every z and k. 
The study of the Schwarzian derivative of / could be a good point of view. Moreover for 
every p £ N, it seems that we can find z and k such that the attracting orbit has 2 p for prime 
period. To conclude, can we write a complete classification of the map with topological 
conjugacy: Two maps f : A ^ A and g : B ^ B are said to be topologically conjugate if 
there exists a homeomorphism h : A ^ B such that f = h~^ o g o h. In this case / and g are 
equivalent in terms of their dynamics. For instance if a; is a fixed point of /, h{x) is a fixed 
point of g. 

6 Appendix 

In this section, we remind elementary definitions and results on discrete dynamical systems we use 
in section [3] and SI For a more thorough presentation, we refer the reader to [3J S] • 

So let M be a metric space and / : M —>• M a mapping. We are interested in the discrete 
dynamical system corresponding to 

f{x) = X, P{x) = o f(x), n £ N, X £ M (6.1) 


Definition 6.1 (fixed and periodic points) A pointy £ M satisfying y = f{y) is called a fixed 
point of f. 

A point y G M satisfying y = f^{y) for p £ N* and y ^ /”(p) for n £ {1, • • • ,p — 1} is called a 
periodic point of f of prime period p. 

Note that a fixed point of / is a periodic point of / of prime period 1. 


Definition 6.2 (attracting and repelling periodic points) Let y G M be a periodic point of 
f of prime period p. 

1 . y is an attracting periodic point of f (of prime period p) if there exists an open neighborhood 
U of y such that 


Vx £ U, lim /"^(x) = y. 

n —>■+00 


2 . 


y is a repelling periodic point of f (of prime period p) if there exists an open neighborhood U 
of y such that 

Vx £ U\{y}, 3n £ N s.t. /"^(x) ^ U. 


Definition 6.3 (hyperbolic periodic point) Lety G M be a periodic point of f of prime period 
p, and let A = D(fP){y) be the Jacobian matrix of fP at y. If A is invertible and has none of its 
eigenvalues on the unit circle, we call y a hyperbolic periodic point of f (of prime period p). 
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Proposition 6.1 Lety be an hyperbolic periodic point of f of prime period p, and let A = D{fP){y) 
be the Jacobian matrix of f^ at y. 

1. If all the eigenvalues of A are inside the unit circle, then y is an attracting periodic point of 
f.yis said to be linearly attracting. 

2. If all the eigenvalues of A are outside the unit circle, then y is a repelling periodic point of 

/• 

Note that a hyperbolic periodic point which is neither attracting nor repelling is called a saddle 
point. In this case, it is of interest to define the stable and unstable sets of the mapping fP. 

Definition 6.4 (stable and unstable sets) Let y be a hyperbolic periodic point of f of prime 
period p. We define the stable and unstable sets W^{y) and W'^{y) as 

( 6 . 2 ) 

(6.3) 

Note that if y is a attracting periodic point of /, then there exists an open neighborhood U of 
y such that W'^{y) (111 = {y} and W^{y) (lU = U. Similarly, if y is repelling, then there exists an 
open neighborhood U oi y such that W‘‘{y) (lU = {y} and lP“(j/) (lU = U. 

Theorem 6.5 (Hartman-Grobman Theorem) Let y be a hyperbolic fixed point of f. Then 
there exists a neighborhood U of y and a homeomorphism h : U ^ M such that 

f\u = h~^ oAoh, 

where A is the Jacobian matrix of f at y. 

This theorem states that in a neighborhood of a hyperbolic fixed point, / is topologically 
conjugated to its linearization. So in this neighborhood, the behavior of the dynamical system is 
qualitatively the same as the one of its linearization. This leads to the stable manifold theorem. 

Theorem 6.6 Let y be a hyperbolic fixed point of f. Then 

1 . W^{y) is a smooth manifold and its tangent space at y is the stable space of the linearization 
of f at y. 

2 . lT“(j/) is a smooth manifold and its tangent space aty is the unstable space of the linearization 
of f at y. 


W^{y) := \x G M \ lim /"^(x) = y 

I n—>-+oo 

lT“(y) := /x € M I lim /"^(x) = y 
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